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I.  INTRODUCTION 


The  electrostatic  potential  due  to  a  planar  triangular  patch  of 
surface  charge  of  unit  charge  density  in  free  space  is  called  V  and  is 
given  by 


V(r ) 


(1) 


Here,  T  denotes  the  surface  of  the  triangular  patch,  ds'  is  a  differential 

element  of  surface  area  on  T,  _r'  is  the  radius  vector  to  the  location  of  ds' 

r  is  the  radius  vector  to  the  point  at  which  V  is  evaluated,  and  e  is  the 
—  o 

permittivity  of  free  space.  The  potential  integral  (1)  arises  in  numerical 
solutions  of  electrostatic,  magnetostatic,  and  electromagnetic  problems  in¬ 
volving  arbitrarily  shaped  surfaces  modeled  by  planar  triangular  subdo¬ 
mains  [1] ,  [2]  . 

The  potential  integral  (1)  has  been  evaluated  by  several  authors 
[1].  [3]-[6],  It  is  perplexing  to  note  that  the  form  of  V  obtained  in  any 
one  of  the  five  references  [1]  and  [3]-[6]  is  different  from  that  obtained 
in  any  of  the  others.  Although  the  form  obtained  in  [6]  is  perhaps  the  most 
convenient  for  numerical  work,  the  objective  of  the  present  report  is  to 
give  a  detailed  derivation  of  the  form  obtained  in  [5].  This  form  is 
[5,  Eq.  (A)]. 


II.  TRANSFORMATION  TO  AREA  COORDINATES 

In  this  section,  the  integral  on  the  right-hand  side  of  (1)  is  ex¬ 
pressed  as  an  iterated  integral  with  respect  to  two  area  coordinates. 


As  shown  in  Fig.  1,  the  triangle  T  has  vertices  r^ ,  r^ *  a°d  £3 • 
In  Fig.  1,  the  point  _r'  on  T  is  called  P',  and  T  is  divided  into  the 
three  subareas  A^ ,  A and  A^ .  To  facilitate  the  evaluation  of  V(r),  area 
coordinates  (£,n.O  are  defined  by 


t. 


(2) 


where  A  is  the  total  area  of  the  triangle  T. 
It  is  now  shown  that 


r'  =  (1-n-O^  +  rl£2  +  C£3 
In  order  to  obtain  (3),  we  express  r'  as 


(3) 


r'=r1+  R'  (4) 

where,  as  shown  in  Fig.  2,  R'  is  the  vector  drawn  from  vertex  1  to  the 
point  P'.  The  dotted  line  in  Fig.  2  passes  through  the  point  P',  is 
parallel  to  the  line  13,  and  intersects  the  line  12  at  the  point  C.  It 
is  obvious  from  Fig.  2  that 


R'  =  1C  +  CP'  (5) 

where  1C  is  the  vector  drawn  from  vertex  1  to  the  point  C  and  CP'  is  the 
vector  drawn  from  C  to  P'.  Equation  (3)  is  recast  as 


R' 


trip 


(6) 


where  1C  is  the  length  of  the  vector  1C  and  CP'  is  the  length  of  the 
vector  CP'.  From  Fig.  2,  it  is  easy  to  see  that 


4 


and 


1C 


V^3 


H1  —3 


CP' 


h 

H 


Hence,  (6)  becomes 


(7) 


(8) 


R'  =  j )  (9) 

Substitution  of  (9)  into  (4)  gives  the  desired  result  (3). 

In  regard  to  the  differential  element  of  area  ds'  in  (1),  Fig.  3 
shows  a  finite  element  of  area  As.  The  area  As  is  the  area  of  the  parallelo¬ 
gram  bounded  by  the  lines  corresponding  to  ",  "+'r,,  r,  and  l+t '  .  The 
vertices  of  this  parallelogram  are  the  points  a(n,  :+AO,  b(n+An»C+A£)  * 
c(n,rj  ,  and  d (r+Ao , A )  •  The  area  of  this  parallelogram  is  given  by 

As  =  | ac  x  ab  :  (10) 

where  ac  and  ab  are  the  vectors  drawn  from  the  point  a  to  the  points 
c  and  b,  respectively.  Specifically, 


ac  =  _tr '  (rvA)  -  r'(p,«;+AA)  (11) 

Substitution  of  (3)  into  (11)  gives 

FF  =  Ar,  (r  j  -  r  3)  (12) 

Similarly , 

ab  =  "r'(r2-r1)  (13) 


Substitution  of  (12)  and  (13)  into  (10)  gives 
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Representation  of  's  as  the  area  of  the  parallelogram 
v.’hose  vertices  arc  a,  b,  c,  and  d. 


3 


Range  of  P  along  a  line  o 


f  constant  £  • 


6 


As  =  |(r2-r1)  x  (r  -r1)  |Ar,An 


(14) 


Thanks  to  (14),  (1)  becomes 


V  (r ) 


X  (— 3  — 1}  1  f 


1  1-C 


4  ft  c 


df, 


dn 

R 


(15) 


where 


0  0 


R  =  ir  -  r’ 


(16) 


and  _r'  is  given  by  (3).  The  limits  of  integration  in  (15)  were  ob¬ 
tained  from  Fig.  4. 

III.  THE  INTEGRAL  WITH  RESPECT  TO  n 

The  integral  with  respect  to  n  in  (15)  is  evaluated  in  this 
section. 

As  the  first  step  in  the  evaluation  of  this  integral,  R  is  expressed 
in  terms  of  n  and  I.  Substitution  of  (3)  for  _r'  in  expression  (16)  for  R 
and  subsequent  squaring  give 


which  becomes 


where 


R2  =  I n (£2”£1 )  +  C(r3-r1)  -  ( r-r  3 ) ! 2 


2  2 

R  =  on  +  2 bn  +  y 


(17) 


(18) 


a  =  |  r  -  | 


f-  =  (r^-r^  •  (^(^-It)  -  (l-.r^))  > 

i 

2  I 


Y  =  ! ',(r3-rr)  -  ( r-r 3 ) 


(19) 


Completing  the  square  in  (18),  we  obtain 


R 


■((';  +  ~-)2  + 


(20) 


Being  the  square  of  the  leneth  of  one  of  the  sides  of  the  triangle 
'■  ’  0.  Now,  thanks  to  (20), 

1- '  1  -  ‘ 


t 

o 


d_n 

R 


d- 


(-  +  V  + 


T, 


(21) 


2 

Since,  from  (17),  R“  ■  0,  it  is  evident  from  (20)  that 


(n  +  :r+  ;■  0 


(22) 


for  all  finite  values  of  '  and  in  particular,  (22)  is  true 

for  T)  =  -?l  ,  so  that 

>,  -  ■  *"  •  0  (23) 

-i 

It  is  evident  from  (19)  that  is  a  quadratic  function  of  q  where 

9 

the  coefficient  of  ^  is 


I  1 2  ,  2  9 

1 12  -  .  |  t;3  -  r  L  j  -  ( (r,  -  r^)  •  (r_3  -  rx) )  “ 

The  above  quantity  is  positive  because  the  vectors  (£0  -  _r  )  and 

(j"3  -  r^),  being  two  sides  of  ttie  triangle  T,  can  not  be  parallel  to 

2 

each  other.  So  far,  it  has  been  established  -ty  -  P  satisfies  the 

inequality  (23)  and  is  a  quadratic  function  of  where  the  coefficient 
2 

of  '  is  positive.  Consequently,  there  can  not  be  more  than  one  value 


of  '  for  which 


8 


<V  -  I'2  =  0  (24) 

Hence,  as  far  as  the  integration  with  respect  to  in  (15)  is  concerned, 
we  may,  when  integrating  witli  respect  to  t;,  assume  that 

-  i-2  -0  (25) 

In  view  of  (25),  application  of  [7,  Formula  200.01.]  to  the 
integral  on  the  right-hand  side  of  (21)  gives 


1- 

f 

o 


d_; 

R 


[  In (*  + 


o  +  V  + 


(26) 


Multiplication  of  the  argument  of  the  logarithm  in  (26)  by  any  posi¬ 
tive  constant  will  not  change  the  value  of  (26).  Choosing  to  multiply 
by  v.i  and  evaluating  t  ho  right-hand  side  of  (26)  at  the  limits  0  and 
1-",  we  obtain 


where 


1- 

=  —J.  [In  (R.,  +  x7)  -  In  (R1  +  Xj)] 

^  V  : 

0 


(27) 


(28) 


Ro  = 


.(1--)  + 


-r 


(29) 


x 


1 


(30) 


x 


(31) 
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According  to  (20) ,  the  argument  of  the  square  root  in  (28)  is 

2  2 

R  at  n  =  0,  and  the  argument  of  the  square  root  in  (29)  is  R  at 

7 

0  =  1-f, .  Now,  use  of  (17)  for  R"  reduces  (28)  and  (29)  to 


where 


r~2 

R  .  =  /  A  .  ■  +  8  .  '  +  C  .  ,  i  =  1 , 2 

i  V  i  i  i 


A. 

l 


(32) 


(33) 


B .  =  -  2 (r  -  r . )  •  (r  -  r  . ) 
x  1  i  i 


(34) 


C  . 

l 


(35) 


Substitution  of  (1°)  for  and  .  in  (10)  and  (31)  gives 


wiie  re 


x .  =  0 . '  +  K  j  ,  i  =  1,2 


II.  =  -  — 


(  r  7  - 


V 


(r3  -  ’V 


-  r 


(r., 


V 


(r 
r . 


r.) 

i 


(38) 


(37) 


(38) 


IV.  INEQUALITIES  FOR  A.,  B.,  C  ,  0.,  AND  E 

_ _ i _ i _ i _  \ _ x 

Because  r^,  r0,  and  r^  are  the  vertices  of  the  triangle  T,  A ^ 
of  (33)  satisfies  the  inequality 

A  0,  i  =  1,2  (39) 

Being  two  sides  of  the  triangle  T,  tiie  vectors  r^  ~  r^  and  r^  -  _r 
can  never  he  parallel  to  each  other.  Hence,  it  is  evident  from  (33) 


and  (37)  that 
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A.  -  D2  0  ,  i  =  1,2  (40) 

For  later  convenience,  we  desire  to  have 

0.  4  0  ,  i  -  1,2  (41) 

i 

Since  the  triangle  T  can  not  have  two  right  angles,  it  is  evident  from 
(37)  that  =  D0  =  0  is  impossible.  If  either  or  is  zero,  then 
we  can  cyclically  permute  the  assignment  of  the  vectors  _r  ^  »  and 

to  the  vertices  of  the  triangle  T  until  neither  nor  is  zero.  Hence 
we  can,  without  loss  of  generality,  assume  that  (41)  holds. 

Tt  is  evident  from  (23),  (28),  (29),  and  (32)  that 

R2  =  A.\2  +  B.'  +  C.  0  ■  4  <  (42) 

ii  i  l  — 

2 

Now,  R^  can  be  expressed  as 

2 

B.  ?  4A.C.  -  B. 

R2  =  A.f('  +  ~-)L  +  — ^ (43) 

1  1  2Ai  4A2 

1 

Clearly,  the  minimum  value  of  R*T  with  respect  to  *,  in  (43)  is  given  by 

4A. C.  -  B2 

min(R7)  =  — ^77 - —  (44) 

l  4A . 

i 

where  min  denotes  minimum  value.  In  view  of  (39),  it  is  apparent 
from  (42)  and  (44)  that 

4A.C.  -  B2  0  (45) 

ii  l  — 

The  definitions  (33)  -  (35)  of  A.,  B  ,  and  C.  reveal  that  equality  is 
attained  in  (45)  only  when  the  point  r  lies  on  the  line  which  passes 
through  the  points  r  and  r\  . 
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Inspection  of  (28)-(31)  reveals  that 

,2 


R.  -  x. 
1  l 


-  _il 


-  B 


(46) 


Taking  the  minimum  with  respect  to  C,  of  both  sides  of  (A6) ,  we  obtain 


„2 


.  /„2  2.  .  ,ny  -  c  . 

min(R.  -  x.)  =  min( - ) 

li  m 


(47) 


Since  a  is  positive,  it  is  apparent  from  (20)  that  the  minimum  of  R  with 


respect  to  q  and  ”,  is  the  minimum  of  (07  -  P  )/<<  with  respect  to  C,. 


,2 


min(R“)  =  min( — - ; — ) 

v  i 

As  a  consequence  of  (47)  and  (48), 


(48) 


2  2  2 

min  (R  .  -  x  . )  =  min  (R  ) 

i  l 


(49) 


Tf  the  point  r  is  not  in  the  plane  of  the  triangle  T,  then 


r  -  tj  can  not  be  expressed  as  a  linear  combination  of  the  vectors 


r2  -  r^  and  r^  -  r ^  which  appear  in  (17).  As  a  result,  the  vector  subject 


to  the  magnitude  operation  in  (17)  can  never  he  zero  so  that 


min(R“)  '  0 


(50) 


As  a  consequence  of  (49)  and  (50), 


2  2 

min  (R  .  -  x . ) 
l  l 


(51) 


The  definitions  (32)  and  (36)  of  R  and  x.  lead  to 


n  0  ,  ,  3.  -  2D.E. 

2  2  2S  .2  ,i  li 

R  .  -  x  .  =  (A .  -D.)(.,  +  - 5 - 

llll 


C.  -  E. 
-i _ b) 

O  f 


(52) 


A.  -  D. 

l  l 


A.  -  D. 

l  l 
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Upon  completion  of  the  square,  (52)  becomes 


B.  -  2D.E  d. 

(A.  -  nf)  [(■  +  -i - +  -i] 

1  1  2 (A .  -  D2)  D2 

11  1 


where 


(B.  -  2n.E  )2 
1 _ 1  1 

2  2 

4 (A .  -  D.) 

l  l 


(53) 


(54) 


In  view  of  (40)  and  (41),  it  is  evident  from  (51)  and  (53)  that 

d2  0  (55) 

The  above  inequality  holds  whenever  the  point  £  is  not  in  the  plane 
of  the  triangle  T. 

If  the  point  £  is  in  the  plane  of  the  triangle  T,  then  r_  -  _r^ 
can  be  expressed  as  a  linear  combination  of  the  vectors  r^  -  and 

r^  -  £^  which  appear  in  (17).  Consequently, 

min(R2)  =  0  (56) 

Retracing  the  development  (50)-(55)  with  (50)  replaced  by  (56),  we 
obtain 

d2  =  0  (57) 

l 

Equation  (57)  holds  whenever  the  point  £  is  in  the  plane  of  the 
triangle  T. 

In  view  of  (19),  substitution  of  (27)  into  (15)  gives 

’  (la  -  I,  >  ■  (£■»  -!,)'?  i  f 

V(r )  =  - y— - r-1 - - - —  l  (-D  [  1  n  (R  .  +  x.)]dr,  (58) 

o  -2  -1  i=l  q 

The  integral  in  (58)  will  he  evaluated  next. 
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Instead  of  the  precise  integral  in  (58),  we  consider,  for 
generality  and  notational  simplicity,  the  indefinite  integral  1^  de- 
f ind  bv 


r 


T1  = 

[  1  n  (V  +  v)  ]  dx 

J 

(59) 

where 

V  = 

\J  Ax2  +  Bx  +  C 

(60) 

V  = 

Dx  +  E 

(61) 

Here,  A,B,C,D,  and  E  represent  the  previously  defined 

constants  A . ,  B . , 

l  l 

C  ,  D  ,  and  E.,  respect ivel v , 
ii  l 

From  (39),  (40),  (41), 

(45),  (54),  (55), 

and  (57),  we  have 

A  0 

(62) 

A  - 

n2  0 

(63) 

D  ^  0 

(64) 

4 AC  - 

B2  ^  0 

(65) 

9 

0 

(66) 

where 

2  2  r  C  - 

d  =  n  [ — — 

E2  (B  -  2DE) 2  i 

2  2  2 

n  4  (a  -  n  ) 

(67) 

A  - 

2 

V.  THE  INTEGRA!.  OF  THF.  LOGARITHM;  4 AC  -  B  =0 

In  this  section,  the  integral  I  of  (59)  is  evaluated  for  the 
case  in  which 

9 

4 AC  -  B“  = 

Equation  (67)  can  be  rewritten  as 


0 


(68) 
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d2  = 


D  "2  2  KA  '  D2)(4AC  -  B2)  -  (2AE  -  BD)2] 


4A(A  -  D  ) 

In  view  of  (66)  and  (69),  it  is  apparent  that  (68)  is  accompanied  by 


(69) 


2AE  -  BD  =  0 


d2  =  0 


(70) 

(71) 


Completing  the  square  on  the  right-hand  side  of  (60),  we  obtain 


v  /,  .  B  s2  ,  4AC  -  B 

V  =  v  A  y  (x  +  — )  +  - ^ - 


(72) 


4  A 


Equation  (68)  reduces  (72)  to 


V  =  lA  *  +  2A 


Substitution  of  (73)  and  (61)  into  (59)  gives 


f  — ,  g  , 

1^  =  [ln(*A  x  +  —  +  Dx  +  E)ldx 


(73) 


(74) 


Assuming  that  x  +  B/(2A)  does  not  change  sign  during  the  integration 
in  (74),  we  obtain  [7,  Formula  610.) 


r  -  (V  +  v)  In  (V  +  v)  -  (V  +  v)  B  n 

1  T  n  2A 

*  A  -  D 


(75) 


t  -  (V  +  v) 1 n (V  +  v)  -  (V  +  v),  .  _B_  n 

,  _  —  2A 

1  .A  +  D  ZA 


(76) 


The  combination  of  (75)  and  (76)  will  be  a  valid  representation 
for  the  integral  (59)  only  if  this  combination  is  continuous  at 


x  +  JK  =  0 
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According  to  (61)  and  (73) 


V  +  .-.A!x+i|+Dh+i)+iM_jJ2 


(77) 


Equation  (70)  reduces  (77)  to 


V  +  v  =  >  A  x  +  +  D(x  +  ^) 


(78) 


It  is  evident  from  (78)  that  V  +  v  approaches  zero  as 


x  +  -*  0 

2A 

Therefore,  the  right-hand  sides  of  both  (75)  and  (76)  approach  zero  as 


*  +  £  -  ° 


Consequently,  the  combination  of  (75)  and  (76)  is  continuous  at 

<  +  H-° 

Hence,  the  combination  of  (75)  and  (76)  is  a  valid  representation  for 
the  integral  (59).  This  combination  is  expressed  as 


2 Ax  +  B  ,  N  B 

I,  =  — — -  ln(V  +  v)  -  x  - 


2  A 


2A 


(79) 


Omitting  the  constant  -B/(2A)  from  (79),  we  obtain 


2Ax  +  B  ,  N 

T  -j  =  ~  2A~"  ln(V  +  v)  _  x 

Expression  (80)  is  1  of  (59)  for  the  case  in  which  (68)  holds. 


(80) 
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VI.  THE  INTEGRAL  OF  THE  LOGARITHM;  4  AC  -  B2  >  Q,  d2  =  0 

In  this  section,  the  integral  1^  of  (59)  is  evaluated  for  the 
case  in  which 

4AC  -  B2  ■  0  (81) 

d2  =  0  (82) 


The  substitution  [8,  Sec.  2.251] 


V  =  t  -  .  A  x  (83) 

is  used. 

Recalling  that  V  is  given  by  (60),  we  square  both  sides  of  (83) 
to  obtain 


Bx  +  C=t^-2iAxt 


(84) 


and  we  find  that 


2.  A  t  +  B 


(85) 


2  (v  At2  +  Bt  +  v  A  C) 
dx  =  - — — ~ -  dt 

(2»  At  +  B) 


(86) 


V  = 


> A  t  +  Bt  +  i  A  C 
2v  A  t  +  B 


(87) 


V  +  v 


G  A  +  D)t  4-  (2i  A  E  +  B)t  +  BE  4-  (>  A  -  D)C 
2*  At  +  B 


(88) 


Upon  completion  of  the  square  in  the  numerator,  (88)  becomes 

,  —  .  r  ,  ,  2/ A  E  +  B  s2  ,  d2(./A  -  D)2, 

G  A  +  D)  [  (t  +  - ^ - )  +  - - - o -  ' 

2  G  A  +  D)  D _ 


2  v  A  t  +  B 


V  +  v 


(89) 
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where 

.2  -  D2  rBE  +  (/A  -  D)C  (2/AE  +  B)21 

d  ~  _  2  _  _  2  (90) 

0  A  -  D)  »  A  +  D  4 (/A  +  D) 

It  is  easy  to  show  that  the  right-hand  side  of  (90)  is  equal  to  that 
of  (69). 

Expressions  (85)  -  (89)  would  not  he  well-behaved  if  the  quantity 
2i  A  t  +  B  could  pass  through  zero.  This  quantity  will  never  be  zero  if 
it  can  be  shown  that 

t  +  JL  ■  n  (91) 

2i  A 

The  inequality  (91)  is  established  in  the  following  manner.  From  (83), 


t  *  V  +  .  A  x  (92) 

Substitution  of  (72)  into  (92)  and  subsequent  addition  of  B/(2>  A)  to 
both  sides  of  (q2)  give 


A  ( 


(93) 


In  view  of  (62)  and  (81),  the  right-hand  side  of  (93)  is  positive. 
Hence,  the  inequality  (ql)  holds.  Consequently,  expressions  (85)-(89) 
are  well-behaved. 

Thanks  to  (82),  substitution  of  (86)  and  (89)  into  (59)  gwes 


where 


fl  =  2 


i  A  t _  + _ Bt  +  i  A  C 

(2.  A  t  +  B)“ 


In  w  dt 


■  A  +  n  +  2,  A  E  +  B^2 

2.  A  t  +  B  2(.  A  +  ID 


(94) 


(95) 
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The  identity 


reduces  (94)  to 


'Ar+Bt  +  >Ai'  1  a r ,>  (96) 


(2,  A  t  +  B)  4i\A 


4A(t  +  — 5— )2 

2/h 


where 


r  =  J-  (i  +  4AC--  1_  x  ) 

1  2.A  2  4A  3 


(97) 


l  =  1 n  w  dt 

2  J 

f  In  w 

S  R  ~2 

3  j  (t  +  -5-r 

2,  A 


dt 


(98) 

(99) 


Use  of  (95)  leads  to 

In  w  =  1  n  (‘  A  *_-)  -  In  (t  +  — -)  +  2  l.n  1 1  +  (100> 

2>  A  2>  A  2(»  A  +  D) 

Thanks  to  the  integration  formula  [7,  Formula  610.1 

r 


In  x  dx  =  x  In. 


J 


the  integral  (98)  becomes 


t  1  n  (*  A  *  ^)  -  w^lnwj  +  +  2w2ln  :W2l 

2  >  A 


-  2w„ 


(101) 


where 


Wj  =  t  + 


W0  =  t  + 


_B 
2\  A 

2t  A  K  +  B 
2  0  A  +  D) 


(102) 


(103) 


Upon  substitution  of  (92)  for  t,  (101)  becomes 


(V+/A  x)ln  (- - — )  -  w^lnw^  +  +  Zw^lnlw^j 

2v'A 


(104) 


where 


2Ax  +  B 


(105) 


_  2Ax  +  B  +  2AE  -  BD 

2TA  2v  "X  (v  A  +  D) 


(106) 


Substituting  (100)  into  (99)  and  changing  the  variable  of  inte¬ 
gration  from  t  to  Wj  of  (102),  we  obtain 


1  ,  A  +  Dn  f  lnWl  ,  ,  „  [  1  ,  ,  ,  2AE  -  BD  | 

I,  »  -  —  In  ( - — )  -  ;  — +  2  !  —  ln  ,wt  +  — ~ — — - 1  dw 


3  w 


1  2>  A 


1  r~  <—  1  1 

2>/A(v  A  +  D) 


The  integrals  in  (107)  are  evaluated  using  the  integration  formulas 
[7,  Formulas  611.2.  and  621.2.1 


1  In  x  . 

■  - —  dx  = 


1  nx  1 


f  lna+bx 


_  -’-1  dx  =  —  Inx  -  (—  H — )ln  a+bx 
2  a  x  a 


The  result  is 


1  ,  ,»  A  +  D,  .  lnwl  ,  1  ,  4>  A  (v  A  +  D)  , 

I  - - In  ( - )  +  -  +  —  +  — =-t-= — -gr — -  In  w 

3  W]  2iJ  w:  wj  2AE  -  BD  1 


_2  ,J_  +  A  0  A  +  D)  ln.  ; 

Z  Cw  2AE  -  BD  ’  1  21 

where  w^  and  w^  are  given  by  (105)  and  (106).  Equations  (69),  (81), 
and  (82)  assure  that  the  denominator  2AE  -  BD  in  (110)  is  not  zero. 
Substitution  of  (104)  and  (110)  into  (97)  gives 


(110) 
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1  r/  ,  4AC  -  B2,  ,  /lr  ,  -  4 AC  -  B2,  ,  ,/a  +  D, 

-  (w  -  2w  +  t — - )  +  (V  +  >  A  x - — - )  In  ( - ) 

2,  A  1  2  "*A  w-[  4A  Wj  ^ 


,  ,  ,  4AC  -  B~  /  1  4v  A  (iA  +  D)U1 

+  (-w  +  - — -  ( —  +  — — l - — — ))ln  w 

I  4A  2AE  -  BD  1 


.  4 AC  -  IT  ,  1  .  2v  A  (»  A  +  D)u  ,  -  .  1 

+  2(W2 - 4A  (^7  2AE-BD— })  ’"‘V3 


(111) 


With  V,  ,  and  given  by  (60),  (105),  and  (106),  simple  algebra  re¬ 
duces  (111)  to 


,  1  ,  .  2 Ax  +  B  2AF.  -  BD  ,  4Ax  +  B  .»•  A  +  I). 

I  =  — - -  [  ( - —— - — — - )  +  — — - —  In  ( - ) 

2*  A  v  A  v  A  (>  A  +  D)  2.  A  2v  A 


,2 Ax  +  B  (4 AC  -  B2)  (.  A  +  D),  , 

-  (— - —  -  - - - - - — - - )  In  w 

>  A  .A  (2AF.  -  BD) 


+  2<“2_±J5  +  — _  OAC  -  l-ljuun),  ln:  I, 

,  A  2.  A  (>  A  +  D)  2v  A  ( 2AE  -  BD)  2 


(112) 


From  (69)  and  (82),  we  obtain 


-  it  -  m 


a  -  n 


Substitution  of  (113)  for  4AC  -  in  (112)  gives 


2»  A  F.  +  B  ,  ,  B  .  ,  ,i  A  +  I). 

I  =  -x - — - - - —  +  ( x  +  — )  In  ( - 

2.  A  (>  A  +  I))  '4;  2  ,  A 


(113) 


2.  A  E  -  B  ,  ,  ,  „  ,  B-2DE  „  ,  , 

(v  -  - - - )  In  w,  +  2(x  +  - - - —  )ln|w2] 


2,  A  (,  A  -  I>) 


2  (A  -  I)  ) 


(114) 


F.xpression  (114)  is  recast  as 
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,  ,  B  -  2DE  w,  /k  +  D,  ,  ,  „  ,  i  u 

I  =  (x  +  - =— )  (In  ( - — — )  -  In  w  +  2  In lw9 | ) 

2  (A  -  D  )  2.  A 

—  2  f- 

2AE  -  BD  ,  2.  A  E  +  B  ,  4ADE  -  AB  -  BD  .  ,v  A  +  D 

+  — — - —  In  w  -  x  - •+ - — = -  In  ( - 

2.  A  (A  -  D‘)  2.  A  (>  A  +  D)  4A  (A  -  DZ)  2v  A 

(115) 

It  is  evident  from  (84),  (82),  (102),  and  (103)  that 


In  (lA_±J2)  _  in  w  +2  ln'w7]  =  ln(V  +  v)  (116) 

2.  A 

Us  ins;  (116)  and  (105)  and  omitting  the  two  constant  terms  that  follow 
-x  in  (115),  we  obtain 


,  B  -  20E  .  s  ,  2AM  -  BD  ,  ,  2 Ax  +  B, 

I  =  (x  +  -- - v)  1  n  (V  +  v)  +  -3 - —  In  (V  +  - — — )  -x 


2 (A  -  D“) 


2i  A  (A  -  D“ ) 


2>  A 


Expression  (117)  is  I  of  (54)  for  the  case  in  which  (81)  and  (82)  hold. 


O 

VII.  THE  INTEGRAL  OF  THE  LOGARITHM;  cT  0 


In  this  section,  the  integral  1  of  (59)  is  evaluated  for  the  case 
in  which 

d2  ■  0  (118) 

It  is  evident  from  (69)  that  (118)  is  accompanied  by 

4 AC  -  B2  0  (119) 

Integrating  (59)  bv  parts,  we  obtain 

I  =  x  1 n  w  +  I ,  (120) 

1  <4 


(117) 


<n*-. 
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where 


and 


w  =  V  +  v 


I, 


2DV  +  2 Ax  4-  B 
2wV 


xdx 


(121) 

(122) 


The  integral  l,  will  he  evaluated  by  means  of  the  substitution 

*4 

(83).  In  this  substitution,  the  new  variable  t  is  given  by 


t  =  V  +  i  A  x 


Expressions  for  x  and  V  in  terms  of  w  and  t  are 


t  -  w  +  E 
.  A  -  1) 


i  A  u  -  Dt  -  . A  E 


»  A  -  D 


023) 


024) 


(125) 


Substitution  of  (124)  and  (1251  into  the  numerator  of  the  integrand 

of  I ,  g  i ve s 
4 


2  (>  A  -  0) 


[2  (.A  +  n)  t  -  2  >  A  W  4-  2,  A  ::  +  Bl(w  -  t  -  e)  d,.  (l26) 
w  V 


which  becomes 


1  , 


*-♦ 


(127) 


where 


f  2(2»  A  +  I))t 


2.  A  w  4-  4>  A  E  +  B 
V 


dx 


(128) 


1 6  = 


ft  A  4-  0 )  t  4-  ( 4>  A  1.  4- 


ini: 

wV 


4-  B)  t  4-  (2v  A  E  +  B)  E 


dx 


(129) 


In  view  of  (61),  substitution  of  (121)  for  w  and  (123)  for  t  in 
(128)  gives 


(130) 
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where 


f 

l  =  2  (v  A  +  D)  dx 


(131) 


I  = 
8 


A  I-.  + 
V 


dx 


(1  32) 


From  (131),  we  obtain 


I  =  2  0  A  +  D)x 


(133) 


Replacement  o'  1?  hv  2B  -  R  in  the  numerator  of  the  integrand 
of  (132)  and  use  of  (60)  vive 


where 


I  +  ( 2*.  A  : ' 


n) 


10 


(134) 


2  +  R>:  +  0 


(1  3  3) 


I 


o 


(1  3b) 


I’erf ormin-.t  the  inti-. -rat  ion  in  (!  !‘>  )  ,  wi  obtain 


0  37) 


Substitution  of  (72)  for  V  in  (!)b)  cives 


dx 


10 


L  'I  -  4A0  -  B 

v/(x  +  ,-■)  + - - 

4A“ 


(1  38) 


Application  of  [7,  Formula  200.0!.]  to  stives 


B  ,  2  .  4 AC  -  IT 


'  1  0 


'  In  (x  +  r.  +,/  (x  +  ...  -  -  )  0  39) 

.  A  “  4A“ 


Multiplication  of  the  argument  of  the  1  ost.ir  i  thn  in  (139)  bv  a  constant 


I 

s 
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is  permissible  because  it  is  equivalent  to  adding  a  constant  to  I-jq- 
Multiplying  the  argument  of  the  logarithm  by  » A,  we  obtain 


r 


10 


(v  + 


2 Ax  + 
2.  A 


Substituting  (137)  and  (140)  into  (134)  and  then  substituting 
(134)  and  (133)  into  (130),  we  obtain 


1  . 
j 


(»  A  +  D)x  +  4V  + 


(2,  A  F. 


ja  (.,  + 


2>  A 


The  integration  in  (129)  will  he  performed  with  respect  to  t. 
After  substitution  of  (86)  for  dx ,  (87)  for  V,  and  (88)  for  w,  (129) 
becomes 

,  ■  2  0  A  +  D)t2  +  (4;  A  K  +  2PK  +  B)t  +  (2>  A  F.  +  B)E  rft 

6  -  (>  A  +  D)t"  +  (2,  A  K  +  B)t  +  BK  +  G  A  -  D)C 

Expression  (142)  reduces  to 

U  =  4t  +  2rn 

where 

t  =  ■  ( 2 OF  -  B)t  +  (2.  A  F.  -  B) F.  -  2  (,  A  -  P)C  _  dt 

11  J  (,  A  +  D)t2  +  (2>  A  E  +  B)t  +  BF.  +  (.A  -  D)C 

Completing,  the  square  in  the  denominator  of  (144),  we  obtain 


1  '  (2DE  -  B)t  +  C  .  A  E  -  B)E  -  2  (>  A  -  D)C 


11 


,  A  +  n  -  (t  +  2.  A  E  +  B}2  +  d20  A  -  D)_ 


dt 


2  (.  a  +  n) 


D 


.2  -  Bv  (no).  Expression  (143)  is  rewritten  as 


where  d  is  given  by 


(140) 

(141) 

(142) 

(143) 

(144) 

(145) 
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‘11 


2DI1  -  B 


a  +  n 


12 


where 


4BP11  +  4C (A  -  T)2)  -  4AE2  -  B2 
2  (>  A  +  D)2 


13 


t  +  2-  A-  A  +  B 


2  (>  A  +  D) 


12 


■sr  dt 


(t  +  - 


9-  A  K  +  B  s  2  .  d^(,  A  -  D)' 


-)  + 


2  (»  A  +  D) 


dt 


13  '  (t  |  -*>  A  E  +  B)2  |  d2(,  A  -  D) 

2  (>  A  +  D)  d" 

From  (90),  we  find  that 


O 

cT  = 


- -  (4Bi)i:  +  4C (A  -  n2)  -  4AF2  -  B2) 

4 (A  -  0  )- 


so  that  (146)  reduces  to 


r  =  2DK  -  B  .  c  "(-  A  -  D)~ 

ii  -  12  2  n 

.  a  +  n  i) 


Recoqniziny  that  the  numerator  of  the  integrand  in  (147)  is 
half  of  the  derivative  of  the  denominator,  we  obtain 


,  -  i  In  [  (t  +  +  O2; 


12  “  2(>  A  +  0) 

liquation  (89)  reduces  (151)  to 


1  ,  r  (V  +  v)(2,  A  t  +  B)  , 

I]2  =  o  n  - J 


n 


v  a  +  n 


Because  t  is  riven  by  (92), 


2.  A  t  +  B  =  2>  A  V  +  2  Ax  +  B 


Equation  (153)  reduces  (152)  to 


(146) 


(147) 


(148) 


(149) 


(150) 


(151) 


(152) 


(153) 
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1  =  \  In  (V  +  v)  +  i  In  (V  +  2AX-+  B)  +  |  In 

2.  A  “  v  A  +  D 

Omitting  the  constant  term  in  (154),  we  obtain 


I 


12 


1  1  >  A'-  +  R 

~  In  (V  +  v)  +  In  (V  +  ~) 

2i  A 


Application  of'  (7,  Formula  120.!.]  to  (148)  cives 

A  F.  -r  B, 


D 


13 


d  (.  A  -  !') 


t  an 


:m  +  -—-■—) 
2  0  A  +  0 )  j 

d(.  A  -  1') 


In  (156)  and  in  the  remainder  of  this  report,  d  is  the  positive  square 
of  the  right-hand  side  of  (  1  )  ( .an!  tin  value  o'  tan  1  is  restricted  t* 
between  -  anti  '/2.  Suhstitul  i  mi  o;  (l)..‘l  lor  t  in  (166)  fives 


n  -i 

=  — —  -tan  u 
<!(.  A  -  !>  i 


where 


( 2_(,  A  +  i))(V_+_LA  x_)  +2w\  I!  +  B)j) 
5  2d  (A  -  ir) 


Thanks  to  (155)  and  (1)7),  (150)  becomes 


I  ,  ■  ■ ■  1 1  n (v  +  v)  +  i„(v  +  SSLil,)  -  “tArOM 

2  (i  A  +  f) )  2>  A 


D 


Substitution  of  (159)  and  (92)  inti'  (143)  gives 


I  =  4V  +  4>  A  x  +  — - - -  [  1  n  (V  +  v)  4-  ln(V  +  -Ax  +  B)] 

6 


i  A  +  n 


2.  A 


id  (i  A  -  D)  .  -1 

- - -  tan  w3 


(154) 

(155) 

(156) 

r  pot 
I  it 

(157) 

(1  58) 

(159) 


(160) 


Subtracting.  (160)  from  (141),  we  obtain 
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I  -  I  =  -20  A  -  D)x  -  — - -  1  n  ( V  +  v) 

5  6  ,  A  +  !» 


2AI-  -  BI>  ,  ,  2 Ax  +  H  ,  4d  (  A  -  D)  _  -1 

— - - -  I  n  (V  + - )  +  — ■  -  - -  tan  w 

•  A  (.A  +  D)  2*  A  ’  3 


(161) 


Thanks  to  (161)  and  (127),  (1  20)  become 


,  ,  ,  B  -  2  HE  ,,  ,  ,  ,  2 -  BI)  ,  .  2 Ax  +  B, 

I  =  (x  +  — - —  )ln(\  +  v)  +  -  — - ,  In  (V  +  - — — ) 

‘  2  (A  -  I'")  2.  A  (A  -  0“)  2*  A 


-  x  +  ,,  tan  Wj 


(162) 


Expression  (162)  is  ot"  (  A))  via-n  o. 


VIII.  AN  F.Ol'  i  VA1.KNT  KXl'K  ESS  1  r>N  !’•  i  HI  INVERSE  TANOENT 


In  this  section,  an  equivalent  expression  is  obtained  for  the 
quantity  tan  h-  that  appears  in  (162).  It  is  evident  from  (157)  that 


,  "I  d(  A  -  0) 

t.,n  W(= - - - 1 13 


(163) 


Backtracking,  we  see  that  I  is  riven  by  (148).  Using  (86),  (87),  and 
(84)  to  chance  (1*8)  back  into  an  integral  with  respect  to  x,  we  obtain 


A  +  0 


dx 


I  i 


V(V  +  v) 


(164) 


In  (164),  V  and  v  are  given  by  (60)  and  (61),  respectively.  Rationalizing 
V  +  v  in  the  denominator  of  the  integrand  of  (164),  we  obtain 


A  +  n 


1  3 


",4  -  l15> 


(165) 


where 


f  _ dx _ 

9  *  9 

-  V4"  -  v~ 


(166) 


A 


Equations  (171)  and  (09)  reduce  (174)  to 


A 


(175) 


Application  of  [7,  Formula  120.1.]  to  (173)  qives 


hs 


n  -l  a  t. 

- - a  Lan  (— H ) 

d (A  -  D") 


076) 


It  is  evident  from  (72)  that 


A  T  =  V 


(177) 


so  that  (1 76)  becomes 


1  5  2 

d  (A  -  D  ) 


tan 


(178) 


Expression  (178)  is  I  of  (167)  for  the  case  in  which  (171)  holds. 

Consider  I , r  of  (167)  for  the  case  in  which 
1  5 


B  -  2 1 )  K  ,  B 

2  f  A 

A  -  I) 


(179) 


In  this  case,  I,-  will  hi'  evaluated  bv  means  of  the  substitution  [8, 

1  :> 


Sec.  2.252] 


St  +  C9 

t  +  1 


(180) 


where  C  and  C.  are  constants  that  will  he  determined  ^iter.  Application 
3  4 

of  this  substitution  to  V  ot  (60)  e.ives 


C3l  +  C4  2  C3t  +  C4 . 

V  =  ■/  A(  +  B (— ~r  -.— -)  +  C 


t  + 


t  +  1 


(181) 


which  becomes 
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V  = 


(AC2.  +  BC_  +  C)t2  +  (2AC  C.  +  BC,  +  BC,  +  2C)t  +  AC2  +  BC,  +  C 
3 _ 3 _ 3  4  3  4  4  4 


(t  +  1)' 

From  (60),  (61),  and  (180),  we  obtain 


V2  -  v2  =  {  [  (A  -  n“)C2  +  (B  -  2DF.)C,  +  C  -  I-2]t2  +  [  2  (A  -  D2)C„C, 

3  3  3  4 

+  (B  -  2DE)  (C  +  C,)  +  2  (C  -  I-:2)  It  +  (A  -  D2)C2  +  (B  -  2DE)C,  +  C  -  E2 
3  a  4  4 


The  constants  C.  anti  C,  are  determined  hv  annihilating,  the 
3  4 

coefficients  of  t  in  (182)  and  (183).  These  coefficients  will  vanish  if 


2AC„  C.  +  B  ( C  +  (',)  +  2C  =  0 
3  4  3  4 


)  I 

ire.,  <:,  +  ni:(c  +  c, )  +  i;“  =  o 

>  )  -4 


Equation  (183)  is  rewritten  as 


(C,n  +  K)  (])<;  ,  +  E)  =  0 

•4 

The  above  equation  will  be  satisfied  it 


Substitution  of 

(187)  for 

in  (184) 

it  i  ves 

BE 

(4  BO 

-  2(10 

-  2  All 

The  inequa 1  i t v 

(179)  assures 

that 

BO  -  2AE  4  0 

With  C  and  riven  bv  (187)  and  (188),  V  of  (182)  becomes 
3  a 


(182) 


7(t+l) 

(183) 

(184) 

(185) 

(186) 

(187) 

(188) 


(189) 
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where 


„  l  /  2  ,  .BE  -  20)72  .  BF.  -  2CjT  7 

=  Tt  +  1  A’L  A  Bl)  -  2AK  +  B(BD  -  2AE)  +  C 


c  =  (cn  -  bde  +  af.2)/d2 


(190) 


(191) 


It  can  be  shown  that 


■  <4AC  -  “",C 


where 


h  =  (BP  -  2AK)/rT 


n“b“ 


(192) 


093) 


Equation  (192)  reduces  (190)  to 


t  + 


, 2  +  (1AC  -  B~)c 
!)“b“ 


liquation  (I  t1))  can  be  recast  as 


094) 


>1'  =  — (4  (A  -  D“V  -  h“V) 

WA  -  i ) '  ) 

It  is  evident  ; re™  (6  3),  lbs),  (lie)  and  (193)  that 


(195) 


e  n 


Thanks  to  (119)  and  (!9b),  V  of  (194)  tan  he  expressed  as 

.  c  fl  2 

4  =  ---f  A  +  q 


where 


d  4 AC  -  B" 

lih 


096) 


097) 


098) 


Sub-.t  iiui  ion  e!  (187)  and  (188)  into  (183)  gives 


I 
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VJ  -  V-  =  •[(A  -  ir)(£)2  -  (B  -  2DE)(£)  +  C-E2]t2  +  (A-I)2)  (|| - 


BD  -  2AE 


BF  -  9(T)  9  ? 

+  (B  -  2DE)(^--^~)  +  C  -  E“ }  /  ( t  +  L)2 


(A  -  !)“>([;■)“  -  (B  -  2I>F)(-|)  +  C  -  F.2  =  c 


099) 


(200) 


(A  -  n2)<f^$r  -  (B  -  +  c-r2  =  (Aac_-b:.)c  -„4c>-  (201) 

■  D  b 

Hence,  (194)  reduces  to 

V"  -  v“  =  -1 — (202) 

(t  +  ir 


where  ~ 

_»  O 
-  ~-- 

and  (|  is  given  by  (198).  Equat  ion  (149)  can  be  recast  as 

•> 

2  1)“  2  2 
d  =  - - — ;r-  (4  AC  -  B  -  41)  c) 

4  (A  -  !)-)“ 

It  is  evident  from  (19K),  (203),  (204)  and  (118)  that 


(203) 


q  p 


From  (180),  we  obtain 


0..  "  C, 

dx  =  — - — iydt 


(t  +  IT 

Substitution  of  (187)  and  (188)  into  (206)  gives 


lib  ( t  +  1  ) 

where  b  and  c  are  given  bv  (143)  and  (191),  respectively.  Substituting 
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(180)  for  x  in  (61)  and  using  (187),  (188),  (191),  and  (193),  we  obtain 


-  2c 

b(t  +  1) 


(208) 


2 

Substitution  of  (197)  for  V,  (202)  for  V 


2 

v  , 


(207)  for  dx. 


and  (208)  for  v  in  (167)  yields 


4i  c  f 

nir 


t  +  i 


-dt 


,  T  2 

t  +  1  (t~  +  q"  -  p 


')  /  9 

“)  V  t~ 


(209) 


+  q 


Assuming  that  t  +  1  does  not  change  sign  during  the  integration  in  (209), 
we  obtain 


where 


>15  =  -  F 


1  5 


(210) 

(211) 


F  =  - 


Ai  C 


( 


dt 


Db2  J  2  2  2. 

(t  +  q  -  p  ) 


(212) 


A2  7 


+  q 


Avoiding  the  question  as  to  whether  or  not  the  combination  of  (210)  and 

(211)  is  a  valid  representat ion  for  I  of  (209),  we  proceed  to  evaluate 

(212)  and  transform  (210)  and  (211)  into  the  domain  of  x. 

The  stibst  i tut  ion 


and  subsequent  equations 


dT  = 


(t- 


2,3/2 

q  ) 


dt 


O  ■> 

l"  +  <f 


(213) 


(21  A) 


(215) 


I 
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reduce  (212)  to 


F  = 


4i  C 


f  dT 


2  2,  2  2 
Db  p  T2  +  3  ~  P 
2 
P 


Application  of  [7,  Formula  120.1.]  to  (216)  gives 


(216) 


F 


-1  /  P'r 
tan  ( — 1 - - — ) 


FI  2 

f  1  -  P 


Db“p  /cf  -  p‘ 

In  (217),  T  is  given  hv  (213). 

With  a  view  toward  expressing  T  of  (213)  in  terms  of  x,  we  use 
(208)  to  obtain 


(217) 


t  + 


bv 


(218) 


It  is  evident  from  (218)  that 


t  =  - 


bv  +  2c 


bv 


Fq nation  (197)  g  i ve s 


f 


t  +  1'V 


(219) 


(220) 


Substituting  (218)  for  t+1  in  (220),  we  obtain 


f 


2  2  2.  c  V 

t  +  q  =  ■  -,--T 
bv  i 


(221) 


Thanks  to  (219)  and  (221),  (213)  becomes 


T  =  - 


bv  (bv  +  2c) 


(222) 


2  bv '  i  c  V 

Before  substituting  (222)  into  (217),  we  use  (198),  (203),  and 
(204 )  to  oh t a i n 

/T~  2  2 (A  -  n2)d 

V'l  -  P  =  — 1 - 

n"  b 


(223) 
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Now,  (222),  (223),  and  (203)  transform  (217)  to 


I’ll  -1  -v  (bv  +  2c)D\ 

- —  tan  ( - - -  -  —  ) 

b  (A  -  lT)d  2 ' v ' dV (A  -  D~) 


(224) 


which  is  equivalent  to 


bvl) 

bv  (A  -  ir’)d 


- ! 

tan  0 


(225) 


where 


(bv  +  2c  )  1) “ 
2dV  (A  -  !)“) 


(226) 


In  view  o!  (19b),  it  is  evident  from  (218)  that  bv  0  when 
t  -1  and  that  bv  0  win  a  t  -1  .  Hence,  with  F  given  bv  (225), 
(210)  becomes 

si  - ; 

I  -  -  - - —  !  ni  1 . ,  bv  •  0 

’  (A  -  0~)d 

and  (211)  becomes 


(227) 


1  =  _ - n- — - —  i  '  c,  «,v  n  (228) 

1  )  _ 

(A  -  n  )d 

The  combination  of  (227)  and  (228)  will  be  a  valid  representation  for 
Ij,_  of  (167)  if  this  combination  is  cont  inuous  at  bv  =  0.  Since  the 
right-hand  sides  of  (227)  and  (228)  are  identical,  the  combination  of 
(227)  and  (228)  is  obviously  continuous  at  bv  =  0.  Hence,  this  combi¬ 
nation  is  a  valid  representation  for  I  of  (167).  This  representation  is 


15 


(A  - 


I)’ )  d 


tan  C 


(229) 


A 
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I  f 


-1 

=  tan 


C 


then 


C,  =  tan 


(230) 

(231) 


Now,  (231)  implies  that 


tan  (--  -  .'.) 


I 

(I 


(232) 


Because  g  lies  between  -  /2  and  '/2,  —  -  g  must  lie  between  0  and 

It  is  now  evident  from  (232)  that 


-I* 

tan 


(  233) 


-1  * 

where  tan  is  the  inverse  tancent  whose  value  1  ies  hetwc  >n  0  and 
From  (233),  wo  obtain 


Substitution  of  (2  30)  lot-  ,•  in  ( :  \ ;  i  -ivis 


(234) 


tui 

Subst  i  tut  in.-  (2  1))  :  or  t  ar. 
-/ 2,  and  usin  -  (226 1,  we  obt.iin 


-  t  (235) 

in  (229),  omitting  the  constant 


1  5 


(a  -  n“)d 


-1 *,2dV(A  -  0"), 
tan  (— - - - 7^) 

(bv  +  2c)  IT 


(236) 


Expression  (236)  is  I  of  (167)  for  t ho  case  ia  which  (179)  holds. 

If  (236)  were  to  reduce  to  (178)  when  (171)  holds,  then  (236)  would 
also  he  T^r  of  (167)  for  the  case  in  which  (171)  holds.  The  following 
reasoning  is  presented  to  show  that  (236)  reduces  to  (178)  when  (171) 


MSSL. 


f 


is  true.  It  is  evident  f rom  (171)  and  (193)  that 


h  =  0 


so  that  (193)  reduces  to 


A  -  1) 


Solving  (238)  for  e,  substitutin'/,  this;  value  of  c  into  (236),  and  using 
(237),  we  obtain  (178).  Therefore,  (23b)  reduces  to  (178)  when  (171) 
holds.  Hence,  expression  (236)  is  t  .  of  (167)  for  both  the  case  in  which 
(171)  holds  and  the  case  in  which  (179)  holds. 

Concernin'/  the  argument  of  tan  '  in  t.169) ,  we  find  that 


(v  +  B - -)  n  =  u 

2  (A  -  i)“ ) 


where 


2  (A  -  I)“) 

where  v  and  b  are  given  bv  (b|)  and  (193),  respectively.  Hence, 
(169)  becomes 


I.,  =  — -  tan  (~) 
U  d  (A  -  ID  d 


Substituting  (236)  and  (241)  into  (163)  and  then  substituting  (165) 
into  (163),  we  find  that,  within  an  additive  constant. 


-I  1  -1  .it,  -1*  ,2dV(A  -  0“),  , 

tan  w_  =  y  [tan  (-)  -  tan  ( - — - =■)  ] 

3  2  d  (bv  +  2c)n~ 
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IX.  A  GLOBAL  FORM  FOR  THE  IN'TEC.RAl .  Of-  THE  _E< >(. AH !  jj  iM 

Kquat  ions  ( !  Oil ,  (2  39),  .uui  (242 )  .lllow  (162)  t.-  be  rewritten  as 


Ij  =  p  in  ln(V  +  v) 


1)0 


- yin  ( V 

A  (A  -  IT) 


2Ax_  + _ 11 , 

2.  A 


IV: 


-1  ,u,  ,  -1  *  2dV (A  -  17).  , 

+  cl  tan  (-)  -  dtan  (— - - - 4)1 

U 


(bv  +  2c )  I)' 


(243) 


Expression  (243)  is  of  (39)  when  <4  0.  In  this  section,  it  is  shown 

that  (243)  is  a  global  form  for  the  integral  1^  of  (59). 

Expression  (243)  will  he  a  global  form  for  the  integral  1^  of 
(59)  if  (243)  is  I ^  of  (39)  in  all  possible  cases.  Only  three  cases  are 
possible,  the  case  in  which  (118)  holds,  the  case  in  which  (81)  and 
(82)  hold,  and  the  case  in  which  (68)  holds.  By  definition,  (243)  is  1^ 
of  (59)  for  the  case  in  which  (118)  holds.  It  remains  to  be  shown  that 
(243)  is  Ij  of  (59)  for  the  case  in  which  (81)  and  (82)  hold  and  the  case 
in  which  (68)  holds. 

Expression  (243)  will  he  I,  of  (59)  for  the  case  in  which  (81) 
and  (82)  hold  if  (243)  reduces  to  (117)  when  (81)  and  (82)  hold.  Setting 
d  =  0  in  (243)  and  not  inn.  that  tan  '  (u/d)  can  not  exceed  72,  we 
obtain 


3 

1  ,  ,  .  bD  ,  .  2Ax  +  B.  „  ,  /n//\ 

1  =  —  [u  ln(Y  +  v)  -  — — - ~~~7  *n  (^  +  - - ~~)  “  Dx  ]  (2a4) 

1  D  2.  A  (A  -  17)  2>  A 

In  view  of  (193)  and  (239),  it  is  apparent  that  (244)  is  equal  to  (117). 
Therefore,  (243)  is  I ^  of  (59)  for  the  case  in  which  (81)  and  (82)  hold. 

Consider  the  argument  V  +  v  of  the  first  logarithm  in  (244).  It 
is  evident  from  (72)  and  (81)  that  Y  0.  Hence,  if  v  ^  0,  then 


V  +  v  0 .  If  v  0,  we  writ e 
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V  +  v  = 


V  -  v 


(245) 


In  view  ot  (82),  comparison  of  ( 1  #5 H )  with  (166)  gives 


V  -  v"  =  (A  -  D“)(x  + 


R  -  2D!;  _  2 

1 

2  (A  -  l)“) 


Equation  (2)4)  rdurcs  (246)  to 


(246) 


V  -  v 


2  (A  -  I)“)u2 


(247) 


Substitution  of  (  2471  into  (245)  tivi's 


V  +  v 


(A  - 

IWV  -  v) 


(248) 


Accord  im;  to  (248),  if  V  +  v  won-  to  approach  zero,  then  u  would 
approach  zero  simultaneously  so  that  the  u  In  (V  +  v)  term  in  (244) 
would  approaeli  zero  harmlessly.  However,  it  u  is  exactly  zero,  we 
must  agree  to  omit  the  u  ln(V  +  v)  term  from  (244)  without  calculat¬ 
ing  1  n  (V  -r  v) . 

Expression  (24  1)  will  he  I  ^  (if  (44)  for  the  ease  in  which  (68) 
holds  if  (244)  reduces  to  (80)  when  (68)  holds.  Equation  (68)  is 
accompanied  hv  (70),  so  that,  according  to  (lC(j), 


h  -  0  (249) 

Thanks  to  (24<t),  t  !i.  ■  i o.  -ft  ic  ietit  of  the  -econd  logarithm  term  in 
(244)  is  zero.  fn for t  mint e 1 v ,  the  argument  of  t  his  logarithm  could 
he  zero.  If.  because  its  onef f  ie ient  is  zero,  we  agree  to  set  this 

logarithm  ter . enl  to  zero  without  calculating  tlic  logarithm,  then 

(244)  reduce-,  t  ,  > 

t  '  (u  !n(V  +  v)  -  Ox)  (250) 
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12Ax+_B)n  =  ii  (241) 

A 

Substitution  of  (241)  into  (240)  yields 

I  =  ln(V  +  v)  -  x  (242  ) 

Which  is  identical  to  (80).  Hone,-.  (34  U  is  I  of  (49)  for  the  case 
in  which  (68 )  holds.  It  is  evident  from  (78)  that,  whenever  V  +  v 

approaches  zero,  2Ax  +  B  also  approaches  zero  so  that  the  logarithm 
term  in  (242)  approaches  zero  harmlessly.  Nevertheless,  if 
2Ax  +  B  is  exactly  zero,  we  must  agree  to  omit  the  logarithm  term 
from  (252)  without  calculating  ln(V  +  v). 

It  has  been  established  that  (241)  is  I ^  of  (59)  in  all  three 
possible  cases,  the  case  in  which  (118)  holds,  the  case  in  which 
(81)  and  (82)  hold,  and  the  case  in  which  (68)  holds.  Therefore, 

(241)  is  a  global  form  for  the  integral  1^  of  (59). 


X.  USE  OF  THIS  I  NT KORA I.  TO  EVALUATE  THE  POTENTIAL 


In  this  section,  the  global  form  (241)  for  the  integral  1^  of 

(59)  is  used  to  evaluate  the  potential  V(r)  of  (58). 

Replacing  the  constants  A,B,0,I>,  and  E  inherent  in  (59)  and  (243) 

bv  A  ,  B  ,  C .  ,  I'.,  and  E.,  we  obtain 
till  i 


1 


1 


(In  (R .  +  x.))d  =  -  [u.  ln(R.  +  x.) 

j  li  D .  l  it 

0 


i 


b  . I) .  2A  +  B. 

— ■-  ' - 1  n  (R  ,  + 1 - -)  -  D. 

l 


2.  A.  (a.  -  1) 

l  l 


2.  A. 


,  u  .  .  2d , R . (a  .  -  1)  1 

,  ,  .  -l.i.  .  .  -  L*  li  l  ,  , 

+  d.tan  (  )  -  d.tan  (— - —  )  ] 

i  d .  l  b .x .  +  2c .  „ 

l  l  l  l  0 


(253) 


where  and  x.  are  given  bv  (12)  and  (16),  respectively.  Also,  d^  is 


the  positive  square  root  of  t  lie  ri'ht-hand  side  of  (54).  Furthermore, 
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a  =  a  /rr 

i  i  i 


h.  =  (B.D  -  2A.i:.)/D7 

i  l  i  ill 


c.  =  (C.D2  -  B.P.K.  +  A.O/n: 

l  li  ill  li  i 


(256) 


,  i 

u .  -  x .  +  — - 

i  i  2  ( a  .  - 


1) 


(257) 


Equations  (255)  -  (257)  are  similar  to  (l^l),  (191),  and  (240), 
respectively.  It  is  evident  from  (195)  that  an  alternative  expres¬ 


sion  for  d .  is 

l 


c.  .  bT 

j-  _  _ i_  _  _ i _ 

i  a .  -  1  .  ,  ,  N  2 

i  4 (a  .  -  1 ) 


(258) 


In  (253),  tan  is  the  inverse  tangent  whose  value  is  restricted  to 

- 1  * 

lie  between  ~"/2  and  II,  and  tan  is  the  inverse  tangent  whose  value 
is  restricted  to  lie  between  0  and  .  If  d  =  0,  then  the  tan  ^  term 
in  (253)  is  to  he  set  equal  to  zero.  If  the  coefficient  of  one  of  the 
logarithm  terms  in  (25!)  is  zero,  < hen  that  logarithm  term  is  to  be 
set  equal  to  zero. 

When  (253)  is  substituted  into  (58),  the  term  in  (253)  cancels 
out  of  tiie  sum  on  i  and  we  are  left  with 


V  (r )  = 


(r0  -  Tj)  •  (r}  -  r 


l5  2  (-1)1 


r„  -  r. 


i=  1 


— —  [u.  In (R .  +  x.)  +  L. 
D .  l  li  l 


1 


_  1  11  •  _  i  *  2d  .  R  .  (a  .  -  1 ) 

+  d.  tan  (.')  -  d.tan  (v~  -1  *■—  r - )] 

i  d  .  i  b ,x .  +  2c .  _ 

i  i  l  l  0 


(259) 


where 


I 
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b  D. 
i  1 


2»  A .  (a.  -  1) 

l  i 


2A . '  +  B. 

In  (R.  +  - - - -) 

1  2>  a  r 


(260) 


XI.  MAN  I PU  LAT  f  ON  OF  THE  POTENTIAL  INTO  THF.  DESIRED  FORM 

In  this  section,  ttie  term  I.  in  expression  (259)  for  the  poten¬ 
tial  V(r)  is  manipulated  into  the  second  logarithm  term  of  [5,  Eq .  (4)]. 
Then,  (259)  is  rewritten  so  as  to  coincide  with  [5,  Eq.  (4)]. 

Expression  (260)  lor  I.,  is  recast  as 


L. 

l 


h.D.  2v  A.  R.  2A.  ■  +  B. 

- 1 — - -  In  (-jr1— 1  +  — ^r^-1)  +  K  . 

2v  A.  (a.  -  1)  Di  Di  11 


L  1 


where 


h.D. 

i  i 


2.  A. 

(7)-T) 
A.  (a.  -  1)  i1 

l  i 


(261) 


(262) 


It  can  he  shown  that 


2A  .  +  B 

— — =  2a.  x.  +  b. 

n .  ii  i 

i 


(263) 


where  x.,  a.,  and  b.  are  given  by  (36),  (254),  and  (255),  respectively 
Substitution  of  (263)  into  (261)  gives 


D  .  2*  A  .  R  . 

1  n  (  (  -  1 - -  +  2a. x.  +  b.))  +  K  (264) 

,  D .  .  D .  ll  l  1 1 


For  the  purpose  of  obtaining  still  another  form  of  of  (260), 
we  call  the  argument  of  the  logarithm  in  (260)  ARC  and  express  it  as 


(R.  + 

l 


2 A  +  B  2 A.'  +  B . 

.  *)  (r  -  -1- - 1  ) 


A. 


ARC.  = 


R  . 


2  A  .  '  +  B  . 

i  i 


(265) 
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With  R  given  by  (32),  (263)  reduces  to 


ARC 


4  A  C 

_ i _ 

4  A  .  ( R  .  - 

i  i 


-  H 
_ i 

2A  -  Vb” 

j _ _ l 

2>  A. 

i 


(266) 


It  4A.C  -  B.  =  0,  then  it  con  be  inferred  from  (70)  and  (193)  that 

b.  =  0  in  which  case  h.  is  to  be  sit  equal  to  zero.  Hence,  it  is 

permissible  to  assume  that  4A.C.  -  P»7  4  0.  in  this  case,  (45)  im- 

111 

plies  that 


4  A  C. 

i  l 


B2  '•  0 


(267) 


Equation  (267)  assures  that  both  the  numerator  and  the  denominator 
in  (266)  are  positive.  Now,  (266)  allows  (260)  to  be  expressed  as 


b.l) 

_ _ i_i 

2.  A.  (a.  -  1) 

i  i 


In  ( 


2.  A  .  R  .  2A  .  ’  +  B  . 

11  1  -1)  +  K„ 


I) . 


D. 


(268) 


where 


K, 


2  i 


h.n.  4A.C.  -  lit 

■  -  -  '  '  In  (  -  -L  - l) 

2.  A.  (a.  -  II  2.  A.  !). 

i  i  ii 


(269) 


Substitution  of  (263)  into  (266)  rives 


I, . 

i 


h.n.  2v  A  .  R 

ii  ,  ,  t  i 

■  -  1:1  (  "n  ’  ' 

2,  A.  (a.  -  1)  i 

i  i 


H  (2a  x. 

ill 

nV 


b.) 

i 


)  +  K„ 


(270) 


It  is  permissible  t o  omit  the  constants  K  and  K ,  .  in  (264)  and 

1  i  2  l 

(270)  because  thev  will  canct  1  out  wlien  the  limits  of  integration  are 

applied.  hiking  (264)  tor  I>.  0,  taking  (27o)  for  1),  0,  and  omitting 

the  constants  K  .  and  ,  .  v.  obtain 
1  i  .i 
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whore  R.,  x^,  b  . ,  c.,  and  u.  arc  given  by  (32),  (36),  and  (254)—  (257) , 
respectively.  Moreover,  d.  is  the  positive  square  root  of  the  right-hand 
side  of  (238).  In  (273),  tan  '  is  the  inverse  tangent  whose  value  is  re- 

—  [  'k 

stricted  to  lie  between  -  /2  .  nd  /2,  and  tan  is  the  inverse  tangent 
whose  value  is  restricted  to  lie  between  0  and  .  If  d.  =  0,  then  the 
tan  ^  term  in  (272)  is  to  he  set  equal  to  zero.  If  the  coefficient 
of  one  of  the  logarithms  in  (272)  is  zero,  then  that  logarithm  term  is 
to  he  set  equal  to  zero.  Except  for  the  facts  that  there  are  no  magnitude 
signs  on  the  argument  of  its  second  logarithm  and  that  its  last  term  con¬ 
tains  an  asterisk,  the  right-hand  side  of  (272)  is  identical  to  the  right- 
hand  side  of  [3,  Kq .  (4)].  There  is  no  need  to  take  the  magnitude  of  the 
argument  of  the  second  logarithm  in  (2/2)  because  this  argument  can  never 
he  negative.  The  second  inverse  tangent  in  [3,  Eq.  (4)]  has  to  be  on  the 
range  from  0  to  in  order  to  pass  continuously  through  '/2  as  b^x^  +  2c^ 


passes  through  zero. 
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